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LAMINAR FLOW ANALYSIS OF FILM BOILING FROM A HORIZONTAL WIRE 
by Kenneth J. Baumelster and Thomas D. Hamill 
Lewis Research Center 


SUMMARY 


The heat -transfer coefficient for pool film boiling from a horizontal wire was 
analyzed by subdividing the wire into many cells and analyzing the behavior of a partic- 
ular cell. A cell consists of a vapor dome and the thin annular film around it, and is 
assumed to be a time -averaged configuration of the actual situation in which vapor domes 
are growing and vapor is escaping. 

If it is assumed that the steady-state equations can be applied to the time -averaged 
configuration, the Navier-Stokes equations with inertia terms neglected, the continuity 
equation, and the energy equation with convective terms neglected can be solved simul- 
taneously to obtain the heat -transfer coefficient for horizontal wires in film boiling. 

Also, the cell size was assumed to adjust itself to maximize the rate of heat transport. 

The heat -transfer coefficient 1^ is given by 


K 


0.373 
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where l is a characteristic length |^g c a/a(p^ - p y )J a the acceleration of gravity, 
k the vapor thermal conductivity, H* an effective heat of vaporization, p^ and p y the 
liquid and vapor densities, respectively, p the absolute vapor viscosity, T w and T gat 
the wire and surface temperatures, respectively, the wire diameter, and g a 
conversion factor in Newton’s law of motion (carried only as a convenience for workers 
in British units - not necessary in cgs system). This expression agrees with the ex- 
perimental data tabulated by Breen and Westwater for wires with 1/ D q ratios of 0. 015 
to 50. 0. 


II 


INTRODUCTION 

Film-boiling heat transfer is important to such diverse fields as cryogenics where 
even a relatively low wall temperature will induce film boiling, or in quenching where the 
high initial temperature of the metal produces film boiling. Film boiling can occur for 
discrete amounts of liquid such as in Leidenfrost film boiling or on submerged heating 
surfaces. In film boiling, a layer of hot vapor supports a layer of dense fluid. Thus, 
in a gravitational field this situation is inherently unstable. For example, figure 1(a), 
displaying a photograph from reference 1, shows the form that the instabilities take. 

At regularly spaced intervals along the wire, the vapor will break through and escape 
under the influence of gravity into the bulk liquid (see ref. 2 for additonal photographs). 

Film boiling, as shown in figure 1, is characterized by a thin annular vapor blanket 
covering the wire, by spatially periodic domes in which the distance between domes is 
nearly constant with time, and by a temporal periodicity - the birth, growth, and 
departure of a vapor dome which represent some average cycle time. This investigation 
was conducted to determine the heat -transfer and geometric characteristics of this 
phenomenon. 

In one of the earliest works on film boiling, Bromley (ref. 3) presented an analytical 
model for film boiling around large tubes. Later, Breen and Westwater (ref. 4) extended 
Bromley's correlation to film boiling from small wires, and concluded that the critical 
wavelength \ Q of hydrodynamic stability theory 


*c = 2 * 
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was an important parameter. (All symbols are defined in appendix A. ) They correlated 
a large amount of experimental data by a semiempirical modification of Bromley’s 
equation: 
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where the modified latent heat H*, which takes into account sensible heating in the vapor 
film, is given by 


H* = H 
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0- 34 C y (T w - T oat ) 


sat' 


H 


fg 


(3) 


In particular, for tube diameters much greater than the critical wavelength, equation (2) 
reduces to 


h oo = 0.59 


" k3p v (P t ~ P v )H * a 
M(T W - T caf )X 
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(4) 


which, except for a constant, is identical in form to Berenson's equation (ref. 5) for a 
flat plate. 

In this report, the laminar film-boiling model is employed. The predicted heat- 
transfer coefficients agree with the tabulated experimental data of Breen and Westwater 
for wire diameters as small as 0. 00053 centimeter and as large as 4. 6 centimeters. 

In addition, in the limiting condition of small wires, the predicted heat-transfer 
coefficient 1^ is proportional to (A c /D 0 ) to the 0. 75 power rather than to the first power 
as used by Breen and Westwater. This agrees with the experimental data and implies 
that, for very small wires, surface tension dominates the gravitational effects. 

The authors would like to thank Professor Westwater for the use of his original data 
plots which had previously been published in reference 4. 


BASIC MODEL AND EQUATIONS 

Figure 1(a) provides justification for the idealized geometric model of film boiling 
shown in figures 1(b) and 2. This model consists of a thin tubular vapor film between 
the wire and the liquid, with vapor escape points indicated by the dome-shaped areas. 
This cellular pattern is distributed in a periodic fashion along the wire with a wavelength 
equal to A. 

The major portion of the heat transport is assumed to occur across the thin vapor 
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film of the tubular portion of the cell. The domes are considered to be so thick that 
essentially no heat is conducted through them. Physically, the domes act as hydro- 
dynamic sinks into which the vapor generated in the annular region is released. A cell 
consists of a vapor dome and an annulus in which heat transfer to the liquid takes place 
by conduction, convection, and radiation through a superheated vapor film; however, in 
this analysis, radiation and convection are neglected. The assumption is made that the 
heat transferred from the wire evaporates the liquid from the interface and that the 
vapor flows in a laminar manner, as indicated by the arrows in figure 2, into the dome- 
shaped sink. 

The vapor dome grows as the vapor enters, breaks away from the wire, and escapes 
into the bulk liquid (see fig. 1(a)). This growth and escape cycle of the vapor dome 
repeats itself. The proposed model assumes the existence of a time-averaged configu- 
ration where all the velocity, pressure, and temperature fields are at steady state. 

Thus, this statistically idealized configuration (fig. 2) represents an average "snapshot" 
of the actual system. 

The following additional assumptions are made in developing the analytical model: 

(1) The inertia terms in the Navier-Stokes equation are negligible compared to the 
viscous terms. This assumption can be justified on the ground that the small thickness 
of the film prevents the buildup of sizeable velocities. Additional support is given to 
this assumption by references 6 and 7, which consider the Leidenfrost film boiling of 
drops. 

(2) In the tubular portion, a uniform vapor gap thickness between the wire and the 
liquid is assumed. Photographic evidence (see fig. 1(a)) justifies this assumption. 

(3) The axial conduction term in the energy equation is assumed to be negligible in 
comparison with the radial term, 


5 2 T « 1 _9_ 
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(4) Sensible heating is accounted for in the final correlation of the heat -transfer 

jjc 

coefficient by replacing the latent heat term with the modified latent heat H that 
was presented by Breen and Westwater (eq. (3)). 

(5) The liquid-vapor interface is at the saturation temperature, and evaporation 
takes place uniformly. 

(6) The velocity and temperature profiles at any instant are assumed to be in steady 
state, and the properties of the flow field are evaluated at the film temperature given by 

T + T . 

T f = — — (6 

1 O 
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(7) The bulk liquid temperature is taken to be at the saturation temperature and 
radiation from the wire to the liquid is assumed to be negligible. 

The governing equations in cylindrical coordinates subject to the preceding 
assumptions are as follows: 


Momentum: 

P v Hr \ Br 2 r 3r r 2 3x 2 j 
- 3P + + 1 3w + 3V\ = 0 

P v \j r 2 r dr 8x 2 J 


Continuity: 


9u 

ar 



= 0 


Energy: 



The boundary conditions for the differential equations are the following: 
At the surface of the wire (fig. 2): 

“N 

r = R 0 

u(R Q ,x) = w(R q ,x) = 0 > 

T(Rq,x) = T w 


(7) 

( 8 ) 


(9) 


( 10 ) 


(ID 


5 


I 


At the liquid interface: 


u(Rpx) = u(R x ) 
w(Rj,x) = 0 
T(Rj,x) = T ga j 


( 12 ) 


From axial symmetry: 


x = 0 


w = 0 from R q r < Rf 


(13) 


In addition, since the values of R^ and u(R^) in the boundary conditions of equation (12) 
are unknowns, an interface energy balance and a pressure balance are also required to 
make the problem determinate. The interface energy balance is given by 


-P v H fg u(Ri) = q=-kiI 


(14) 


r=Ri 


and the pressure balance by 


/ a P(R 1 ,x)27tR 1 dx = 27rR 1 x a P £ 


(15) 


where P_ is assumed to be the average pressure in the vapor at the liquid-vapor 
interface resulting from the weight of the supported liquid, atmospheric pressure, and 
surface tension effects (see appendix C). Equation (14) expresses an interface energy 
balance which assumes that all the heat conducted through the vapor vaporizes the 
liquid. 

In addition, since the two geometric parameters A and x o are unknowns two 

a 

additional assumptions are needed. 
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The first assumption 


lim Jl = 1 (16) 

Rq - * -00 h °° 


expresses in mathematical form the experimental observation of Breen and Westwater, 
as discussed in the introduction, that the heat -transfer coefficient for large wires can be 
calculated by equation (4). As is shown herein, this assumption is, in effect, equivalent 
to using experimental data to determine an unknown constant in the analysis. The 
second assumption 
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represents a hypothesis put forward by the authors in reference 8 that the cell wave- 
length X adjusts itself to maximize the heat removal rate from the surface. In film 
boiling, the liquid-vapor interface is a free boundary that can be deformed in such a way 
as to optimize (in this case, maximize) the heat removal rate from the wire. 

From a mathematical point of view, the problem is now completely formulated. For 
each dependent variable (unknown), there is an appropriate equation or constraint. From 
this point on, the solution of the equations is merely a matter of integration and some 
algebraic simplifications. The details of the analytical procedure are given in appen- 
dixes B to F. However, the general method and principal results of the analysis are 
discussed briefly in the next section. 


DISCUSSION OF ANALYSIS 

The flow equations, momentum and continuity, were converted into a single fourth- 

order partial differential equation (eq. (B3)) by the introduction of the stream function 

and this equation was then solved by separation of variables. A velocity distribution 

o 

was found in terms of one unknown constant of integration £ and a dimensionless 
parameter A, where 
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The constant of integration /3 2 was evaluated by a static pressure balance giving an 
evaporation velocity of the form 


u(A) = 
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The energy transport in the vapor film must now be considered. In terms of the 
dimensionless coordinate r] given by 


R 0 

The energy equation (eq. (10)) becomes 



Equation (20) has the solution for the temperature 


(19) 


( 20 ) 


T = T 
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and for the temperature gradient 


9T _ T w ~ T sat 
dr\ r] In A 
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Combining the expression for the temperature gradient with the expression for 
u(A) in the energy balance at the liquid-vapor interface (see appendix D) gives 


16 k(T. 
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This equation relates the parameter A to the fluid properties and cell geometry. The 
parameters X and P are unknowns, whereas the liquid and vapor properties, as well 

a 

as the radius of the wire, are, of course, known quantities. 

Attention is now directed to obtaining an expression for the heat -transfer coefficient 
from the wire in terms of the quantity A. In appendix E, the heat -transfer coefficient 
(eq. (E8) is shown to be 


»w = 


Rq In A 


(E8) 


where 1^ is the overall heat -transfer coefficient from the wire and 'J> (Rq, £ a ) repre- 
sents the fraction of the total area in which the heat is being conducted into the liquid. 
The first term in equation (E8) expresses the heat-transfer coefficient for a purely 
conducting hollow cylinder. The fraction term ^(Rq, | a ) is added to indicate that the 
vapor dome decreases the effective heat -transfer area. 

Maximizing the heat-transfer coefficient with respect to X according to equa- 
tion (17) (see appendix E) now yields a value for the cell wevelength X of the form 


X 

2 


(1 - { a )aRo 



where the characteristic length l is defined as 


l = 
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The model proposed for this analysis is strictly applicable to boiling from small 
wires. As shown in appendix C, the model does not apply for length to diameter ratios 
less than 0. 35. However, the functional form of the analytical result was such that for 
length to diameter ratios less than 0. 35 the data could be conveniently correlated by 
setting one product of ^ and £ a equal to a experimently determined constant. 

HEAT -TRANSFER COEFFICIENT 

The following two approximations are convenient to use in determining a heat- 
transfer coefficient because of the complications arising from the £ factor: 

cl 

(1) The expression on the left side of equation (D4) is approximated by the expression 

(A 2 - l) 2 + (1 - A 4 ) In A = - - (A - l) 4 (23) 

3 

This expression results from a Taylor Series expansion about A = 1 in which terms 
beyond the fourth power were neglected. As shown in figure 3, this approximation is 
reasonably accurate over a range of A from 1 to 2.4. This range represents those 
values usually seen in practical experiments. 

(2) The expression for the heat -transfer coefficient given by equation (E8) is 
written in approximate form by letting 

In A =* (A - 1) (24) 


which leads to 




k 

R 0 (A - 1) 




(25) 


where 6 is the vapor-gap thickness. 

When these approximations (eqs. (23) and (24), as well as eq. (16)) are used, the 
expression for the heat -transfer coefficient simplifies to the following form (see appen- 
dix F for details) 




= 0. 373 
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T(A,R Q ,l) = 


3^6 
21 A 


z 2 a 2 AR 0 A; 




6A 2 Rq 
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(26) 


(27) 


The latent heat of vaporization H^g was replaced by the modified latent heat of vapori- 
zation H* (eq. (3)) to take into account sensible heating effects. 

In general, the unknown value of A in the preceding equations is a function of the 
fluid and system properties and must be evaluated from the governing equations. How- 
ever, since £„ was not evaluated, some approximate value of A based on experimental 

d. 

data must be chosen. Convenience dictates that one average value of A be chosen over 
the range of practical interest. The value of A becomes significant in the expression 
for the heat-transfer coefficient only for small wire diameters. For large wires, the 
terms containing A in equation (26) become small relative to 1 and the value of A, 
itself, approaches 1. However, for small wires, this is not the case. For example, 
the value of A estimated from figure 1(a) is near 2. The choice of an average value of 
A to be used in equation (F15) was made by a comparison of theory with the experimen- 
tal results. 

The heat-transfer correlation (eq. (F15)) was plotted with Breen and Westwater's 
tabulated data in figure 4(a) for various values of A. An arbitrary choice of A = 1, the 
lowest physically acceptable value of A, best correlates the data. Thus, the equation 
for the heat-transfer coefficient is 


h w = 0. 373 


k H*a (p z - p y )p^ 
^ T w - T sat)* 


1 1/4 


pl/4 
1 1 


(28) 


where 
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As shown in figure 4(b), the correlation in equation (28) fits the data reasonably well. 
The correlations of Bromely, and Breen and Westwater are also shown. The data in 
figure 4 are for 10 liquids (isopropanol, Freon 113, benzene, carbon tetrachloride, 
ethanol, nitrogen, pentane, oxygen, water, and helium), as listed in reference 4. 


SIMPLIFIED EXPRESSION FOR HEAT-TRANSFER COEFFICIENT 

The expression for the heat -transfer coefficient (eq. (28)), simplifies when r is 
expressed in terms of its limiting values for small and large wires. For large wires, 
equation (30) reduces approximately to 

(31) 

L l 


while for small wires 



(32) 


where the subscripts L and S indicate large and small, respectively. Substituting 
these expressions into equation (26) and setting A equal to 1, as before, give for large 
wires 
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1 + 


3 l 
2^6 Rq 


whereas for small wires 


r s 4± t J_/i.V 

b yfi Rq 3^ \Ro/ 


Talcing the value of T as the sum of its asymtotic expressions gives 


r A^ 
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1 + _9__L + J_(_L\ 3 
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(33) 
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(35) 


Substituting the asymtotic value of T given by equation (35) into the expression for 
the heat -transfer coefficient, equation (28), yields 
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= 0. 373 


k H a(Pj - p y )P A 


M(T, 
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Figure 4(c) compares the exact expression for the heat -transfer coefficient, equa- 
tion (28), with the approximate expression, equation (36). Both equations correlate the 
data; however, equation (36) is the simpler to use. 


EFFECT OF ACCELERATION ON FILM BOILING 

Multiplying by g/g each term in equation (36) where the acceleration a appears 
and rearranging give 


r,.3 TT * 


h 7 = 0. 373 
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k°H*(P z - P y )P v g 


M( t w T sa P ^g 
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f— V 
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where 




s c a 


“ 1 1/2 


f <Pl ~ P V >J 


(38) 


and g is a reference acceleration of gravity. 

The acceleration effect on the heat -transfer coefficient is closely related to the wire 
diameter. For small length to diameter ratios (large wires), the heat-transfer coeffi- 
cient varies nearly as the 0. 375 power of a; while for medium length to diameter ratios, 
where the second term in equation (37) dominates, the heat-transfer coefficient varies 
as the 1/4 power of a. For the limiting case of very thin wires, where the third term 
in equation (37) dominates, the heat -transfer coefficient is independent of the accelera- 
tion a. Equation (37) indicates that no conclusions about gravitational effects on film 
boiling are possible unless the wire diameter is also considered. 

Care must be taken when considering ratios of a/g greater than 1. For the case 
of a small wire where the length to diameter ratio is large, the third term in equation 
(37) dominates, and the heat -transfer coefficient is nearly independent of gravity. How- 
ever, a sufficient increase in a would increase the effect of the first two terms in 
equation (37). 

Figure 5 shows the heat -transfer coefficient for various values of a/g on a 0. 476- 
centimeter tube in a Freon 113 film -boiling system for the data of Pomerantz (ref. 9). 

As shown in figure 5, the theory agrees closely with experiment, except at low values 
of a/g. However, the experimental data could conceivably be slightly high for low 
values of a/g, since the heat -transfer coefficient appears to converge to a value of 1. 1 
instead of the expected value of 1. Correlations by Pomerantz (ref. 9), Breen and 
Westwater (ref. 4), and Bromely (ref. 3) are also shown in figure 5. 


LARGE WIRE STABILITY 

In appendix C, the analytical model was shown to be limited to wire diameters such 
that 1/ D q > 0. 35. This limit is very interesting in light of the observation by Breen 
and Westwater (ref. 4) that the heat -transfer coefficients for film boiling are not 
monotonic functions of the tube diameter. Figure 6 shows a plot of the heat -transfer 
coefficient as a function of tube diameter, as given in reference 4. A possible explana- 
tion for the dip and rise in the heat -transfer coefficient is a transition from annular flow 
associated with thin-wire film boiling to circumferential flow associated with large-wire 
film boiling. Thus, in the vicinity of 1/ D 0 equal to 0. 35, the heat -transfer coefficient 
could rise very slightly for increasing wire diameters because of a shift in the flow 


14 


pattern, as illustrated in figure 4(e). The circumferential flow is probably turbulent in 
nature . 


CONCLUSIONS 

A semitheoretical expression for the heat-transfer coefficient associated with pool 
film boiling from a horizontal wire was derived with the assumption that the system 
adjusts itself to maximize the heat transfer. The derived expression agrees with all 
available experimental data, for wire diameters as small as 0. 00053 centimeter and as 
large as 4.6 centimeters, as shown in figure 4(c). For very small wires, the theoreti- 
cal expression for the heat -transfer coefficient increases to the 0. 75 power of length to 
diameter ratio. In addition, the heat -transfer coefficient is shown to be a function of 
gravitational acceleration to the 0. 375 power for large wires and to the zero power for 
very small wires. The gravitational relation of the heat -transfer coefficient agrees 
with the experimental data of Pomerantz for large diameter wires. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 24, 1967, 

129-01-11-02-22. 
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APPENDIX A 


SYMBOLS 


A arm heating area of annular cell 

region, cm^ 

A ce n total heating area of cell, 

27tRqA, cm 2 

a acceleration of gravity, 

cm/sec 2 

a, . arbitrary constants 

Cj, . . . , 5 arbitrary constants 

C y vapor specific heat at con- 

stant volume, cal/(g)(°K) 

D q wire diameter, cm 

d distance, cm (see fig. 8) 

2 

E differential operator defined 

by eq. (B4) 

F(|) arbitrary function 

f arbitrary function (see 

eq. (B5)) 

G arbitrary function (see 

eq. (B5)) 

g reference acceleration of 

gravity, cm/sec 2 

g c conversion factor in 

Newton’s law of motion, 
1.0 g-cm/(sec )(dyne) 

(or 32. 17 (lb mass) (ft)/ 

(lb force) (sec 2 ) for 
British system) 

H* modified latent heat of 

vaporization defined by 
eq. (3), cal/g 


Hfa) arbitrary function 

Hjg latent heat of vaporization, cal/g 

h heat -transfer coefficient, 

cal/(cm 2 )(sec)(°K) 

h exp experimentally determined heat- 
transfer coefficient, 
cal/(cm 2 )(sec)(°K) 

heat -transfer coefficient from 
wire, cal/(cm 2 )(sec)(°K) 

h^ heat-transfer coefficient from 

large wire defined by eq. (4), 
cal/(cm 2 ) (sec) (°K) 

k vapor thermal conductivity, 

cal/(cm)(sec)(°K) 

l characteristic length given by 

eq. (El 3), cm 

l characteristic length evaluated at 

a = g (see eq. (38)), cm 

O 

P absolute pressure, dyne/cm 

P B pressure beneath wire at liquid - 

vapor interface, dyne/cm 2 

P 7 liquid pressure at liquid-vapor 

6 2 
interface, dyne/cm 

P pressure at apex of vapor dome, 
dyne/ cm 

P g average system pressure in vapor, 

dyne/ cm 2 

P 1 pressure at base of vapor dome, 
dyne/ cm 
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^ann 

rate of heat transfer through 
annular region, cal/sec 

r i 

dimensionless function defined 
by eq. (29), r(A = 1) 

Q cell 

rate of heat transfer through 
annular region, cal/sec 

A 

ratio of liquid -vapor -interface 
radius to wire radius , Rj/Rq 

q 

heat flux at liquid-vapor inter- 

6 

vapor gap thickness, R^ - Rq, cm 


face, cal/(sec)(cm^) 

V 

dimensionless radial coordinate, 

R b 

bottom radius of vapor dome, cm 


r/Rfl 


(see fig. 7) 

A 

function defined by eq. (27), cm*^ 

R eq 

equivalent radius of vapor dome, 
cm (see fig. 7) 

X 

X c 

cell wavelength, cm 
critical cell wavelength defined 

R o 

radius of wire, cm 

t/ 

by eq. (1), cm 

R 1 

radius of liquid -vapor inter- 
face, cm 

M 

absolute vapor viscosity, 
g/(cm)(sec) 

r 

radial coordinate, cm 

V 

kinematic viscosity of vapor. 

T 

temperature, °K 


cm 2 /sec 

T f 

film temperature, °K 


dimensionless axial coordinate, 

T sat 

liquid saturation temperature, °K 


x/(A/2) 

T w 

u 

w 

X 

x a 

wire temperature, °K 
radial velocity, cm/sec 
axial velocity, cm/sec 
axial coordinate, cm 
half length of annular region in 

^a 

p l 

p v 

a 

dimensionless annular length, 
x a /(A/2) 

Q 

liquid density, g/cm 

O 

vapor density, g/cm 
surface tension, dyne/cm 

/3 2 

cell, cm 

constant of integration given by 

* 

area correction factor defined by 
eq. (E8) 

3 


2 

eq. (B45), cm /sec 

* 

stream function, cm /sec 

r 

dimensionless function defined 

Superscript: 


by eq. (F8) 

0) 

derivative with respect to inde- 

r A 

asymptotic approximation for 
^i, eq- (35) 


pendent variable 
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APPENDIX B 


SOLUTION OF FLOW EQUATIONS 
Momentum Equations 

The momentum equations (eqs. (7) and (8)) can be transformed into one equation 
with one unknown by introducing the stream function i//(r, x) in which the continuity 
equation (eq. (9)) is satisfied identically and the equations of motion (eqs. (7) and (8)) 
combine to eliminate the pressure terms 


r 3x 


(Bl) 


and 


1 d\p 

w = J- 

r 3r 


(B2) 


thereby leading to a single fourth-order equation of the form 


E V = 0 


(B3) 


where 




9r 2 r3r 3x 2 


(B4) 


If the velocity profiles are assumed to be separable, 


P = f(r)G(x) 


(B5) 


equation (B3) takes the form 


f" M G _ 2 f'"G + — f"G - — f’G + 2f"G M - 2 f T G M + fG"" = 0 

r 


r r 2 r 3 


(B6) 
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Now, a set of functions f and g must be chosen which will satisfy the partial 
differential equation (eq. (B6)) and the boundary conditions. The boundary conditions 
that apply to the momentum equations (eqs. (11) to (13)), can be expressed in stream 
function notation as 


or in terms of f and G 


= R 0 


3x 3r 


(B7) 


r = R. = Const Mi = 0 

1 3x 3r 


x = 0 


= 0 for Rn < r < R- 
3r -u _ _ i 


(B8) 

(B9) 


r = R q fG’ = f'G = 0 


r = R x fG' = Const f’G = 0 


x = 0 f'G = 0 for Rq < r < Rj 


(BIO) 

(Bll) 

(B12) 


The functional forms of f and G cannot be obtained directly from equation (B6) 
because the functions cannot be separated. However, physical intuition and similarity 
techniques used in references 6 to 8 suggest a functional form in which the radial compo- 
nent of velocity is independent of the axial position 

G = x (B13) 

This functional form is consistent with the assumption of a constant vapor gap between 
the liquid and the wire and the assumed boundary condition that the vaporization velocity 
u(R^) is independent of axial position. 

Using equation (B13) reduces equation (B6) to the form 


£T»tf 




(B14) 


Multiplying through by r^ results in 
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r 4 f"" - 21%" + 3r 2 f" - 3rf’ = 0 


(B15) 


which is a fourth-order Cauchy equation (or equidimensional linear equation) having a 
solution of the form 

2 4 2 

f = aj + a 2 r + a^r + a^r In r (B16) 

Therefore, the stream function defined by equation (B5) takes the form 

2 4 2 

4/ = + a 2 r + a g r + a 4 r In r) (B17) 

The analysis was simplified by introducing the following dimensionless variables 

I =— (B 18 ) 

X 

2 

f] = — (B19) 

R 0 

Therefore, the stream function becomes 

P=^-(C 1 + C 2 77 2 + C 3 ? 7 4 + C 4 77 2 In 77 ) (B20) 

2 


Velocity Distribution 

Substituting equation (B20) into equations (Bl) and (B2) results in velocity distribu- 
tions of the form 


1 /Ci 3 \ 

u = — I — + C 2?7 + C 3 rj° + In 17 1 

R 0 \ v ) 


(B21) 
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and 


X 

w = — [2C 2 + 4C 3 r? 2 + C 4 (2 In v + 1)] (B22) 

t ->2 


The boundary conditions given by equations (11) to (13) become 


-3 

1! 

M- 

c 

il 

o 

w = 0 T = T w 

(B23) 

rj = A u = u(A) 

w = 0 T = T sat 

(B24) 

£ = 0 w = 

0 for 1 < 77 < A 

(B25) 


where 


A = 



(B26) 


Applying the boundary conditions, equations (B23) and (B24), to equations (B21) and (B22) 
results in the following matrix equation: 





(B27) 


Solving this matrix equation by Gauss elimination allows the four constants to be ex- 
pressed in terms of one constant of integration, called /3 2 . 


0 2 (ln A + 1 - A 2 ) 

i - 

(B28) 

1 In A 


_ /3 2 (A 2 - 1 - 2 In A) 

(B29) 

In A 
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c 3 = ^ 


c . 2^(1 - a 2 ) 

4 In A 


(B30) 

(B31) 


Therefore, the velocity distributions given by equations (B21) and (B22) become 


u = 


-Z3 2 


In A + 1 - A 2 ( (A 2 - 1 - 2 In A)„ , ^3 , , (1 _ A 2^ In n 

In A In A 


RqL V l 11 A 


(B32) 


and 


- 


X _ 


w = 


Rr 


2 (A 2 - 1 - 2 In A) + ^2 + 2(21n t? + 1)(1 - A 2 ) 


In A 


In A 


(B33) 


The constant 


P 


is determined in the section Static Pressure Balance. 


Pressure Distribution 


The following forms of the momentum equations result from substituting equations 
(B18) and (B19) into equations (7) and (8): 


3P 

dr] 


M 9 2 u 1 9u 

^**0 3 ? j 2 V 



(B34) 


3P 

94 



a?w + iaw + M 

^2 r] 8rj lx I ^2 
\2 / 


(B35) 


Substituting the values of u 
(B34) and (B35) gives 


and w from equations (B32) and (B33) into equations 
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ap - M P 2 L 4(1 - A 2 ) 

9*7 g c RQ Rq L r? In A 


0P 


-16|i(^) 2 /s 2 ! 
8 C R 0 


Integrating equations (B36) and (B37) yields 


p = JLE. 

ScR 2 


4rj 2 + 4(1 - A 2 ) *2^2- 
In A 


+ F(£) + C 5 


—8/a/ — \ 


P = 


‘(l) 




+ H(tj) + C 5 


Therefore, equating equations (B38) and (B39) results in 


P = 




g c R o 


4 r? 2 + 4(1 - A 2 ) 

ln A 



The constant Cg was determined from a simple pressure consideration. The 
pressure boundary condition given by equation (C3) is 


P(A, { ) = P o + + (1 - { a ) ip v A- 

<1 - U- 2 g c 


Evaluating equation (B40) at 77 = A and | = £„ gives 

a 


P(A,U = J^_ 

d n 

e c R S 


2-i 


4 - 8 qf' 

K o 


+ a 


(B36) 

(B37) 

(B38) 

(B39) 

(B40) 

(C3) 

(B41) 
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Combining equations (C3) and (B41) and solving for C 5 yield 


C 5 = P o + 


2 a 


+ (1 


a 


< l - «a> 7 


fa) 2 Pv g 



(B42) 


Static Pressure Balance 


The theoretical velocity and pressure distributions derived earlier specify these 
fields with one unknown, /3 2 . Consequently, an additional mathematical constraint is 
necessary to determine the velocity and pressure fields. 

Substituting equations (B18), (B19), and (B26) into the static pressure balance 
(eq. (15)) yields 

P(A, £)d£ = £ a P g (B43) 


Substituting the pressure distribution, equation (B40), into equation (B43) and collect- 
ing terms give 


0 


2 4/x 


s c R 0 



(B44) 


Substituting the value of C g from equation (B42) and equation (C6), which relates P g 
and P Q , into equation (B44) give 


0 


2 


d-yf -p,>- 


2a 


(1 - U A AR o 



(B45) 


The velocity and pressure distribution are now determined explicitly in terms of the cell 
geometry parameters X and £ a - 
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Interface Velocity 

This analysis predicts a unique relation between the rate of evaporation and the 
thickness of the vapor gap. This relation arises by substituting 77 = A into equa- 
tion (B32) 


u(A) = 


/3 2 (A,! a , A) 
RqA In A 


[(In A)(l - A 4 ) + (A 2 - l) 2 ] 


(B46) 


o 

where the constant /3 is given by equation (B45). Thus, to support the liquid above the 
wire, an injection velocity u(A) at the liquid-vapor interface is required. Combining 
equations (B45) and (B46) gives 




2 O' 


u(A) = -L 
R n 


(i - Q - 
a 2 


AR r 


16p 

3 E C R 0 


V2R r 


.2 

5 a 


(In A)(l - A 4 ) + (A 2 - l) 2 


A In A 


(B47) 
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APPENDIX C 


ESTIMATE OF PRESSURE AT BASE OF VAPOR DOME 

The geometric model shown in figure 2 was redrawn in figure 7. where certain 
pressures are denoted at key points on the model. The vapor dome was assumed to be 
spherical as a first approximation. This approximation is fairly accurate for small 
wires, but not for large wires since the flow tends to be in the vertical direction rather 
than in the axial direction as assumed in this model. However, the problem resulting 
from the change in flow direction was overcome by application of equation (16). 

If P Q is the pressure in the liquid at the apex of the vapor dome, the pressure at 
the base of the vapor dome, shown as point Pj in figure 7, is given by 

P 1- P 0+ 2<T < cl > 
(i - U - 2 g c 

a 2 

Both the pressure jump across a curved interface and the increase in pressure due to a 
vapor head were considered in equation (Cl). 

If it is assumed that the pressure drop in the vapor dome is very small (second 
order), 


P(A, | a ) = P x (C2) 

Therefore, 

P(A, 5 a ) = P 0 ♦ - 2,7 — ♦ (1 - « a > £ P v - 5 - (C3) 

(1 - « a > * 2 g c 

a 2 

This expression is used to determine one constant of integration in the pressure distribu- 
tion function of appendix B, equation (B40). 

Similarly, the pressure in the liquid at the annular liquid-vapor interface is given 
by 


p *= p o + < 1 -«a>f'>*f 


(C4) 
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If the pressure jump across the curved liquid-vapor interface is taken into account, the 

average pressure P in the vapor at the annular liquid-vapor interface is defined as 

s 

p.-p 0 p *f- + ;r (c5) 

2 g c R 1 


The pressure in the annular vapor region is higher at £ = 0 (centerline of symmetry) 
than at £ = £ a (base of vapor dome). This pressure difference represents the driving 
potential for the flow. Therefore, pressure continuity across the interface to some 
fixed reference requires that, in reality, an annular region of constant radius cannot 
exist. However, the assumption of a constant annular vapor region (R^ = constant) 
makes the problem tractable. Consequently, the pressure P Q represents an average 
vapor pressure between £ = 0 and £ = • 

cL 

Expressing Rj in terms of A and Rq by equation (B26) yields 




(C6) 


This expression is used in the static pressure balance, equation (B43), to determine 

o 

the integration constant 0 (eq. (B45)). 

Although the pressure P g and P(A, £ a ) were obtained at the upper surface of the 
annular vapor region, the pressure distribution throughout the annular region at any 
value of £ may be considered constant as a first-order approximation, thereby limiting 
the solution to relatively small wires. The average pressure difference between the top 
and bottom of the annular region is given as 


p s - p b . 2P vRl f- (CD 

g C 

that is, the pressure difference between the top and bottom interfaces can only differ by 
the vapor heat existing between the two extremes. Thus, for small wires with large 
liquid heads 


P b a P s (C8) 

Clearly, the use of P g in the pressure balance (eq. (B43)) should be accurate at 
least to a first-order approximation. For large wires, the failure of equation (C8) is of 
no concern, since the geometric model used in figure 2 breaks down and equation (16) is 


27 



I 


Mill HUB I 


II II 


■■ I I 


1 1 


■ I I 


introduced to overcome this difficulty. 

Some additional interesting facts concerning the shape of the liquid-vapor interface 
of film boiling on a small wire can be obtained by considering the pressure between the 
top and bottom liquid-vapor interfaces. In figure 8, the dome region of the geometric 
model is again shown but in greater detail. For the purposes of the following analysis, 
the dome is considered to have a shape with an equivalent radius of curvature R at 
the top apex and an equivalent radius of curvature at the bottom apex. The assump- 
tion that the top dome is hemispherical is not required here. 

The pressure at the point P B1 calculated along path 1, which goes through the 
vapor, equals 


P Bl = P o 

eq g c B c 


(C9) 


while the pressure at point Pgg calculated along path 2, which goes through the liquid 
equals 


p Bi! = p o tR e<frf ♦ dp J7- + | Z <C10) 

g c s c 

However, from the requirement that the pressure at a given point be independent of the 
path to that point, the ratio of the radius of curvature at the bottom of the dome R^ to 
that at the top of the dome can be estimated. 

Taking the points B1 and B2 to be at one fixed position, the continuity of pressure 
requires that 


P B1 “ P B2 


Substituting equations (C9) and (CIO) into equation (Cll) yields 


Rv. 


R 


eq 


V R eq + d)(p l - p v> ~ 


1 - 




2a 


(Cll) 


(C12) 


Therefore, the radius of curvature beneath the wire is always greater than the 
radius of curvature at the top of the wire, that is, 
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(Cl 3) 


R b >R eq 

Figure 1(a) shows this in an actual physical situation. This result, of course, applies 
only to stationary vapor bubbles, since bubbles rising through a liquid are subject to 
dynamic flow effects. 

In addition, the ratio given by equation (C12) must remain positive for physical 
systems. Thus, 


R eq< R eq + d)(p * " P v> 


2 a 


< i 


Assuming that d is to be small compared with requires that 




^g r 


\Pl ~ P v )a . 


1/2 


(C14) 


(Cl 5) 


If, the radius of the wire is such that 


R 0 > 1.41 


°g r 


ll/2 


Pi ~ Py)* J 


(C16) 


multiple bubbles might form along the wire since R will be less than the diameter of 
the wire. Thus, the flow would begin to shift from the axial to radial type. Therefore, 
an estimate for the lower limit of axial flow is given by 


— >0.35 
D„ — 


(C17) 


where l is defined by equation (El 3). 
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APPENDIX D 


COUPLING OF MOMENTUM AND ENERGY RELATIONS 


The momentum and energy relations are coupled by equation (14) which is repeated 
here in terms of the dimensionless coordinate rj: 


- p y H fg u(A) = 


k 3T 



(Dl) 


Substituting the value of the temperature gradient from equation (22) and the value of 
u(A) from equation (B46) into equation (Dl) gives 


p y H fg /3 2 (A, | a , A) 
RqA In A 


[(A 2 - l) 2 + (1 



k < T w - T sat> 
RqA In A 


(D2) 


Rearranging equation (D2) yields 

(A 2 - l) 2 + (1 - A 4 )ln A = - k(Tw : V . 

P v H fg ^ 2 (X,^ a ,A) 


(D3) 


The complete expression for the interface energy balance was obtained by substituting 

o 

the value of p (eq. (B45)) into equation (D3), which yields 


16 k(T w - T sat )„(-A-r (f 


V2R 


0' 


(A 2 - l) 2 + (1 - A 4 ) In A = 




§c 


2ct a 


(1 - U - AR 0 
a 2 


(D4) 
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APPENDIX E 


MAXIMIZATION OF ENERGY TRANSPORT 

A relation for cell wavelength A., which until now has remained an unknown parame- 
ter, was obtained by using the maximization postulate, equation (17). First, however, 
an expression for the heat -transfer coefficient was needed. 

The total energy transported from the wire to the liquid in a particular cell is 
given by 


^ann ”^ann 3r 


=-kA an 1 H 
r=R 0 R 0 *1 


(El) 


T7=l 


where the annular heat -transfer area is given by 


A ann A cell ^ R 0’ ^a) 


(E2) 


with 


A ce ll = 2ffR 0 X (E3) 

The function ^(Rq, £ a ), represents that fraction of the total area of the cell in which 
heat is being transferred by conduction. This function can be considered as a correction 
factor applicable to the large wire problem in which multiple vapor domes appear. 

Substituting the expression for the temperature gradient at the wire surfaces 
(eq. (22)) into equation (El) results in 


Q 


k(T w - T ( 


sat^ann 


ann 


Rq In A 


(E4) 


where Q ann is the local energy transport based on the annular area adjacent to the 
vapor dome. However, from an experimental point of view the heat -transfer rate and 
the heat -transfer coefficient are based on the area of the wire. Therefore, the heat 
transfer in terms of the total cell area is defined by 

^cell = *\v( T w " T sat^ A cell ( E5 ) 
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but the calculated Q. 


heat -transport rates are taken to be equal, that is, 


ann 


^cell ^ann 

Substituting equations (E4) and (E5) into equation (E6) and rearranging terms yield 


(E6) 




kA 


ann 


R 0 A cell ^ A 


(E7) 


or, using equation (E2), yield 


K = 


Wgp. S a > 

Rq In A 


(E8) 


Thus, equation (E8) is an egression for the heat -transfer coefficient for which the 
maximization hypothesis, equation (17), can be applied. 

Substituting equation (E8) into equation (17) and differentiating with respect to X 
yield 


ax 


(E9) 


Differentiating equation (D3) with respect to X and using equation (E9) yield 


ax 


(E10) 


Performing the preceding differentiation on /Sr (eq. (B45)) and rearranging terms yield 


-d -U-(p? -pJ- + 


6a 


2a 


S, 


c (1 - P) * AR 0 


fi + JL^ = 0 

V 2A dX/ 


(Ell) 


According to equation (E9) dA/dX is zero; thus, solving for the cell wavelength 
yields 

2 "1 
6a2r 2 

1 + — 1 - 1 1 (El 2) 


2 (1 - « a )AHo [_\ l 
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l = 




a (Pi 


^\P_i 


1/2 


(El 3) 


o n 

Since 9 ri/BX' 4 is less than zero, the wavelength defined by equation (El2) is the cell 
radius that maximizes the heat transfer across the liquid-vapor interface. This optimum 
arises because heat transfer by conduction directly under the dome is very small. Thus, 
the dome radius cannot be very large. On the other hand, if X is very large, a large 
pressure difference is needed to force the generated vapor into the dome. However, a 
large pressure difference requires a large dome radius (see appendix C). Thus, it is 
more efficient, from the point of view of maximizing the heat transport, to have many 
cells than to have one massive cell (see fig. 1(a)). There is also a limit, however, on 
the smallness of a cell, since, for infinitely small cells, the surface tension becomes 
important. If the cell wavelength X approaches zero, the pressure rise across the 
vapor dome 2a/R e ^ would become infinite. The pressure in the dome would be greater 
than that in the surrounding liquid; consequently, the dome could not function as a 
hydrodynamic sink. 
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APPENDIX F 


HEAT -TRANSFER COEFFICIENT 

A simplified expression for the heat -transfer coefficient was derived by using two 
assumptions (eqs. (23) and (24)). Substituting equation (23) into equation (D3) and 
rearranging terms give 


(A - 1) = 


3k < T w - T sat) 


L4p v H fg /3 2 (M a , A)J 


1/4 


(FI) 


This expression represents an interface energy balance between heat conduction and 
heat loss due to evaporation. 

Substituting equation (FI) into equation (25) gives the coefficient for heat transfer 
from the wire as 


K = *( R o’Q\ 


4k 3 p y H*/3 2 
_ 3 < T w ' T saP R 0j 


1/4 


(F2) 


where H fg is replaced by H*, given by equation (3) to take into account sensible effects, 
as discussed in the introduction. 

Some mathematical manipulation is now required to determine the final form of h^. 
For convenience, A is defined as 


A = 





Thus, substituting equation (F3) into equation (E12) yields 


X 

2 


(1 - S a > 


A 


(F3) 

l 


(F4) 


In addition, equation (B45) can be written as 
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0 2 = 


3a (Pi ~ P v ) R q 


16 ill 


(i - y - 1 

a 2 




1 - 


21 ‘ 


V 


M a ) 2 | AR od-y 


(F5) 


Substituting equation (F4) into equation (F5) yields 


/s 2 = 


3a (Pi - P v )Rp 
16pZ 


_1_ 

ZA 


[i — ? — + — L 


z 2 a 2 


AR„ 


'1 -L 


=>a 


(F6) 


For mathematical convenience, F(6) is rewritten as 


/3 2 = 


rw 1 - s A r i 

fa(P; - P v )r£] 

^ fi _ 2 + 1 

L 48 V ^a /_ 

1 

3. 

1 

2ZA \ z 2 a 2 ARqA/ 


(F7) 


Defining 


2ZA V Z 2 A 2 AR 0 A z 


(F8) 


gives the expression for as 


0 2 = 



a (Pi ~ P V ) R 0 

L 48 V *a / _ 

MZ 


(F9) 


Substituting this expression for p into the equation for the heat -transfer coefficient 
(eq. (F2)) yields 


K = *<&0’Q 


Ik A - *a\ 2 ] 

1/4 

k 3 a (Pi -P v )H*p y 

_ 36 V / _ 


_ w < T w - T sat>_ 


1/4 


r l/4 


(F10) 


In applying equation (16), the experimentally verified value of h^ of Breen and 
Westwater (eq. (4)) was used. In terms of Z, equation (4) becomes 
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h =0.373 

OO 


k 3 H*a(p z - p y )p v 
. ^ T w - T sat^ . 


1/4 


(Fll) 


Because 


lim A = 

Ho-” ' 


and 


lim T = 1 

R o-°° 


it is necessary that 


£ a ) 


Vi Z 1 - 


,2-i 


36 V £ 


a 


1/4 


= 0. 373 


Assuming that this condition holds for all values of Rq gives 


k 3 H*a {p t - p y )P A 
,^ T w ' T sat^ 




= 0. 373 


,1/4 


,1/4 


(FI 2) 


(FI 3) 


(F14) 


(F15) 
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Figure 4. - Correlations for film boiling on horizontal wires. Wire diameter, 0.00053 to 4.6 centimeters. (Data from ref. 4.) Each type of data point represents a set of test 
conditions within the following ranges: Fluids, water, helium, oxygen, pentane, nitrogen, ethanol, benzene, Freon 113, isopropanol, or carbon tetrachloride; temper- 
ature, 72° to 2600° F; heat-transfer rate 17 to 2120 Btu per hour per square foot per °F ; wire diameter, 0.00022 to 1.895 inches; critical wavelength, 0.068 to 0.660 inches; 
length-diameter ratio, 0.135 to 50. 
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Figure 5. - Influence of acceleration on film boiling in Freon 113 from a horizontal 
0.476-centimeter tube (ref. 9). 
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Figure 6. - Heat -transfer coefficients for film boiling from a 
horizontal tube (ref. 4). 
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Figure 7. - Pressure points in flow model. 



Figure 8. - Calculation paths for determining pressure in vapor dome. 
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